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Abstract 

An innovative application of existing technology is pro- 
posed for attenuating the effects of transient phenom- 
ena, such as rotor-stator and rotor- strut interactions, 
linked to noise and fatigue failure in turbomachinery 
environments. A computational study was designed to 
assess the potential of passive porosity technology as a 
mechanism for alleviating interaction effects by reduc- 
ing the unsteady lift developed on a stator airfoil subject 
to wake impingement. The study involved a typical high 
bypass fan stator airfoil (solid baseline and several 
porous configurations), immersed in a free field and 
exposed to the effects of a transversely moving wake. It 
was found that, for the airfoil under consideration, the 
magnitude of the unsteady lift could be reduced more 
than 18^^ without incurring signiticant performance 
losses. 

Nomenclature 

A Co'oss sectional area 
LH Leading edge 
M Local Mach number 
Re Reynolds number 
Sn Sirouhal number, Sn = ((od)A 27:11) 

T Local temperature 
U Lreestream velocity 
V Wake translational speed 
c Airfoil chord 
c j A i r fo i 1 d rag c( >e ft ic ien t 

C| Airfoil lift coeflicient 
Cp Airfoil pressure coeflicient 
d Diameter of wake generator 
p Local Mow static pressure 

s vSolidity. seconds 

t time 

u Local Mow velocity normal to porous surface 
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X Axial or chord wise distance 
(j) Contraction coeflicient 
y Ratio of specilic heals, y = 1.4 for air 
p Local Mow density 

Subscripts 

1 Upstream of porous surface 

2 Minimum area of Mow through porous surface 

4 Lully mixed Mow region, dowmstream of surface 
c Closed (area) 
max Maximum value 
min Minimum value 
0 Open (area) 
t Stagnation condition 
Freesiream quantity 

0 Start of wake traverse, I = 0.0 seconds 
s Solid stator 

w Wake effects 

Introduction 

Porous media, such as screens and perforated plates, 
have been widely used in Muid dynamics studies since 
the 1940s'. These applications mostly occurred in the 
areas of production /reduction of turbulence and the cre- 
alion/elimination of laige scale velocity or pressure non- 
uniformities. Important examples are the u.se of screens 
to minimize Mow turbulence in wind lunnels% applica- 
tions in boundary layer controL \ and reduction of aero- 
dynamic forces. In the latter application, the porous 
medium is in the form of a perforated outer surface, 
exposed to the Mow at any given incidence. The system 
also includes a solid inner surface such that the volume 
between the two forms a plenum that is lilled with the 
same Muid Mowing over the outer surfaceL This arrange- 
ment, known as passive porosity, redistributes pressure 
on the outer surface by establishing communication 
between regions of high and low pressure through the 
plenum. The pressure redistributi(m, which is associated 
with a minute transfer of mass into and out of the ple- 
num, changes the effective aerodynamic shape of the 
outer surface. Passive porosity has been used to reduce 
wav^' drag and delay buffet onset in transonic airfoils by 
reducing shock-boundary layer interaction'’ L and to 
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reduce large side forces in bodies of revolution at high 
angles of attack caused by asymmeu ic pressure loading 
due to How separation*^. Passive porosity has also been 
used to reduce blade-vortex interaction (BVl) in 
rotorcraft*^. In turbomachinery, it can be applied to the 
reduction of phenomena associated with unsteady (low 
such as rotor-stator and rotor-strut interactions, which 
aie directly related to fan noise and can cause high cycle 
fatigue failure. 

The nature of the unsteady (low in turbomachines is usu- 
ally determined by circulation effects, blade thickness 
effects, and wake effects’-*'^’. The circulation effect is 
caused by the presence of circulation around the moving 
wake-generating blades, and it can be neglected when 
compared to wake effects if the distance between the 
wake generator and the stator is more than 107r - liYA 
of the stator chord. Blade thickness effects are caused by 
the interaction between potential (inviscid) How fields 
around rotor blades of linite thickness. This effect is 
important when the distance between the rotor and sta- 
tor is less than approximately one airfoil chord. The 
wake ( viscous) effect is caused by wake-stator interac- 
tions that form as the wakes shed from the upstream 
rotor impinge on, and then convect past, the downstream 
stator. Since the wake rate of decay is much slower than 
that of the potential How interactions, wake effects can 
still be signiticant several chord lengths downstream. 
Note that this effect would be zero in an inviscid Iluid, 
since wakes are direct consequences of the boundary 
layer development on the generating body. These three 
effects contribute to momentary changes in local How 
incidence and magnitude that give rise to variations in 
surface pressure (i.e., lift). However, the wake effect is 
by far the most important t)f the three. 

The aerodynamic parameters associated with rotor-sta- 
tor and rotor-strut interactions can be grouped into those 
that ( I ) would reduce the wakes shed by the rotors, and 
(2 ) would reduce the response of the stators and struts to 
impinging wakes. Any approach designed to alter these 
groups must do so within the coniines of acceptable per- 
formance losses. For turbomachinery applications, the 
oveiTiding concern is the minimization of flow separa- 
tion on rotor and stator surfaces. I^assive porosity tech- 
nology can be used to modify both groups, the first by 
reducing wake velocity defect and width through reduc- 
tions in rotor drag, and the second by reducing the sta- 
tor/strut response through reductions in fluctuating lift. 
The latter application is the subject of the present study. 


Porous Surface Model 

The porosity model used in this investigation was devel- 
oped in reference 12 from normal screen loss analytical 
and experimental results presented in reference 13. If 
the perforated plate is considered to produce a series of 
jets that coalesce gradually in the downstream direction, 
three uniform states can be determined (Figure 1): the 
free stream well ahead of the screen ( I ). the location of 
minimum area at which the jets formed by the screen 
holes are fully contracted but essentially undiffused (2), 
and the free stream well beyond the screen where the 
flow is fully mixed (3). The upstream and downstream 
areas ai'e equal. The ai*ea at 2 is reduced by the solidity s 
and a contraction coeffleient ()). Thus, 

A, = (1) 

At = Aj())( I s) (2) 

where 


The contraction coeflicient <|) is given by* * 


(t> = + ij (4) 

where the empirical incompressible component 0o is 



I hole 

centerline 


iMgurc 1 .- MtxJel of flow through a single scieen hr)ie. 

The governing equations are derived from conservation 
of mass and momentum laws for the steady, one-dimen 
sional, isen tropic (adiabatic and reversible) How of a 
perfect gas. It is also assumed that there is no total pres- 
sure loss during jet formation. This assumption is valid 
as long as the fluid passing through the hole has an 
inviscid core. Thus, the flow normal to the porous sur- 
face is given by 

Conservation of mass: 

(pu)j = (pu)t())( 1 -s) = (pu)^ (6> 
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Conservation of momentum: 


Pi( 1 + yMf) = 

ii = p,( 1 +YM 5 ) 

(7) 

p^l 1 4- yM 5 ({)( 1 ” s 1 


Adiabatic How; 


^ 1 1 " ^ [2 * 

(8) 

Zero total pressure loss: 


Pil “ Pi2 

(9) 


The model is implemented as a surlace boundary condi- 
tion with the aid oiTwo additional assumptions: (irst, the 
process from state I to state 3 takes place over an intini- 
lesimal distance; second, the pressure inside the plenum 
is constant. As a result, only the How side of the porous 
surface needs to be gridded. and all the information nec- 
essai’y to describe the process is obtained from the 
known mass How in the computational domain and a 
given plenum pressure. 

Assuming that a pressure gradient exists across the 
porous surface, the local How will enter the plenum it 
Pdo„.,,„ > Ppionun.. and exit the plenum il' the converse is 
true. The plenum pressure at every global iteration is 
found by averaging the known surface pressures in the 
computatitinal domain over the area covered by a given 
porous patch. As the solution converges, the plenum 
pressure reaches a steady value for which the net mass 
tlux across the porous surface approaches zero. Detailed 
descriptions of the model implementation and validation 
are given in reference 14 

Computational Tools 

blow Solver 

CH.3D'^ developed at the NASA Langley Research 
Center, was used to calculate the results presented in this 
paper. The code solves the three-dimensional, time- 
dependent, thin-layer approximation to the Reynolds- 
Averaged Navier-Slokes (RANS) equations using a 
Unite vc)lume formulation in generalized coordinates. It 
uses upwind-biased spatial ditlerencing with Roe s tlux- 
difference splitting (FDS) method’^’ for the inviscid 
terms, and central differences for the viscous and heat 
transfer terms. The code, which is second-order accurate 
in space, is advanced in lime with an implicit three-fac- 
tor approximate factorization (AF) scheme. Temporal 
subiierations with muliigrid are used to recover time 
accuracy lost as a result ol the AF apprt)ach during 
unsteady calculations. The pseudo-time subiteration (i- 
TS) method is used in the present work. In this option, a 
pseudo- time term is added to the time-accurate Navier- 
Stokes equations. This extra term is sub-iterated until 


acceptable convergence is obtained, (iood convergence 
was reached using live subiteraiions tor each physical 
time step. C-FI.3I) includes several grid connection stiat- 
egies. a vast aiTay of zero-, one-, and two-equation tur- 
bulence models (linear as well as nonlinear), numeious 
boundary conditions, and translating grid capabilities 
suitable for turbomachinery applications. The results 
presented here were obtained using the two-equation k- 
0) SST model of Menter' \ 

Structured Grid 

The problem of trying to numerically simulate the aero- 
dynamic environment within a turbolan is inherently 
diflicult. However, it can be greatly simplified by con- 
siderine a system composed ol a single rotor blade (or 
any other suitable wake generator) transversely moving 
past a single stator vane. Wake cttecls would still be the 
main cause of unsteadiness within such a system. Thick- 
ness effects, although still present, would be small in 
compaiison to wake effects il the distance between the 
rotor and the stator leading edge is approximately 40^;^ 
of the stator chord Circulaticm effects would be insig- 
nificant at this distance. Furthermore, three-dimensional 
effects such as blade tip and cross-llow ellects can be 
avoided if the system is made two-dimensional. Thus, 
the computational setup chosen for the study consists ol 
a single stator airfoil immersed in a tree (low field, sub- 
ject to the effects of a transversely moving wake. Sche- 
matic representations ol the grid, airloil. and wake are 
given in Figure 2. 

The 9-/one. 2-D grid contains approximately 2.39x10^ 
points/plane, extending to LSc above and below the air- 
foil, 0.26c upstream and 2.5c downstream. Far held 
boundary conditions are applied at the external taces ot 
the grid, interface boundary conditions are applied at all 
the internal faces, and viscous and porous surlace 
boundary conditions are applied at the airloil surlace. 
The vertical extension of the grid, in conjunction with a 
stretched distribution of points toward the tar field, are 
required in order to dissipate retleciions from the bound- 
ai'ies. The majority of the grid points are contained in 
zones 1, 4, and 3. as is necessary to ensure adequate 
propagation of perturbation waves. 

The chosen airfoil was obtained Irom surlace measure- 
ments taken at the mid-span of one of the stator vanes 
composing the Pratt and Whitney Advanced Ducted Pro- 
peller high power fan model currently at NA.SA I.aRC*’^. 
Leading and trailing edges were not defined. Thus, they 
were created by carefully titting circular arcs to the 
eiven discrete points. .Since the section would no longer 
function as part of a stator vane, its incidence had to be 
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adjusted in order to preclude (low separation over its 
surface. Results from an angle of attack sweep indicated 
that there exists a narrow incidence range within which 
the tlow would remain attached over the entire airfoil 
surface. The mid-point of that range. was selected. 


siator 




sumpliiii! 
plane . 


average wake 


f igure 2.- Computational grid. 

Because a strong wake was desired, a cylinder (not 
shown) with a diameter of approximately 0.05c, and 
located 0.4 1 c upstream of the stator lAi plane, was used 
initially as a wake generator. Although a satisfactory 
wake was obtained, its effect on stator lift was heavily 
influenced by the vorticity resulting from the Karman 
vortex street chai*acteristic of How over blunt bodies. 
The effect of the shed vortices manifested itself as oscil- 
lations in C| with frequencies corresponding to Sn - 0.2. 
In order to eliminate these unwanted oscillations, an 
average protile was obtained trom several temporal sam- 
ples of the wake generated by the rod as it approached 
and passed the stator. This average wake, which has a 
?>y/f velocity defect at the plane of the stator LE, repre- 
sented tremendous savings in computing time, since the 
/ones detining the rod would no longer be necessary. 


The average wake, now a user input, is delincd at the 
sampling plane (lace jl ot zone 5, 0.15c downstream of 
the original rod) using a velocity distribution boundai*y 
condition. Through the use of dynamic patched interfac- 
ing, this wake is displaced a distance of approximately 
4c, from face k I to face kmax of zone 5. by a downward 
translation of /ones 5. 8 and 9. At the end of the speci- 
fied traverse, the algorithm translates the moving blocks 
and their solution back to their original position. 

To minimize numerical dissipation as the wake traverses 
the stator flow field, the axial distance between the plane 
of the wake and the stator F.H is kept as small as possi- 
ble, and the cell size ratio between the zones delining 
the wake and the stator does not exceed 3: 1 during the 
entire wake traverse. A translational speed (V) of 78c/s, 
con*esponding to an angular velocity of 5()() rpm and a 
mid-span radius of 1.5c, was chosen for the wake. The 
entire wake traverse requires 2232 time steps. The given 
translational speed. u*aversed distance, and free stream 
Mach number (0.166) guarantee that, at any given time, 
the stator is subject to the passage of a single wake only. 


Results and Discussion 

Hffect of Moving Wake on Airfoil Surface Pressure 

The effect of the moving wake on the llowlield around 
the stator leading edge is sequentially depicted in Fig- 
ures 3a through 3d. CoiTcsponding surface pressure 
coelTicient distributions ai'e presented in Figure 4. In 
general, the pressure on the convex surface of a stator 
vane (or rotor blade) is relatively low, and the pressure 
on the concave side is relatively high. For this reason, 
the convex and concave surfaces are usually called the 
suction and pressure sides, respectively, of the stator. 
The process described here is common in turbomachin- 
ery 

At the start ot the traverse, i = ().() seconds (Figure 3a), 
the wake is located approximately 2c above the stator. 
This relative placement precludes any signilicant inter- 
action between the two. The associated surface pressure 
distribution is as would be expected for the airfoil 
immersed in a free stream at the chosen incidence. Note 
trom Figure 4 that surface pressure is minimum at the 
40V( chord location, coincident with the point of maxi- 
mum airtoil thickness. The recompression of the surface 
pressure downstream of the maximum thickness loca- 
tion indicates that the How remains attached over the 
entire airtoil surface. Close examination of the tlow past 
the Th revealed that the ’‘blip" in the Cp distributions 
results from a combination of airfoil incidence and TF 
geometry, and is not an indication of How separation. 
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The poinl ol' maximum pressure is located at the stator 
KH. Because of the airfoil positive incidence, a point ol 
local minimum pressure, a direct consequence ol local 
How expansion, is located on the suction (leewai'd) side 
at x/c - O.OO.r 

The How around the stator remains virtually unaffected 
by the moving wake until approximately t = 0.0176 s 
(Figure 3b). After this time, the close proximity of the 
wake to the stator LH causes a gradual decrease in the 
local How incidence. At time t = 0.0229 s, the wake 
intersects the stator LH (Figure 3c) causing a large 
reduction in ACp that extends over 90^;^ of the airtoil. 
The change in relative wind incidence shifts the stagna- 


lion point from the origin (x/c = //c = 0) to the suction 
side of the airfoil, and the poinl of local minimum pres- 
sure from the suction side to the pressure (windward) 
side. In the immediate vicinity of the wake, as it passed 
over the airfoil, the velocity was seen to decrease and 
then increase on the suction side, and, conversely, to 
increase and then decrease on the pressure side. As a 
consequence, the pressure rises and then tails on the 
suction side, and falls and then rises on the pressure 
side. This behavior is clearly seen in the surtace Cp dis- 
tributions of Figure 4. The process is then reverted as the 
wake moves past the stator airfoil (Figure 3d). 





l-iaurc 3.- Mach Number contours for solid stator = 0. 166. Re = 1.1 25x 10^\ a - 6^. V - 7Xc/scc. 

Medium grid level. 
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I’igurc Concluded. 



Figure 4.- Surface pressure coenicient distribution for solid 
stator aiiibil. medium grid level. = 0. 166, Re = 1 . 125x 10^’. 


HtTcet of Moving Wake on Airfoil I .ift 

The unsteady nature of the lift generated by a solid sta- 
tor when subjected to wake convection is depicted in 
Figure 5. Note from the rigure that the lift coefricienl 
decreases as the wake approaches the stator, reaches a 
minimum at the point of closest approach between the 
wake and the stator TH (t = 0,0229 s), and then recovers 
as the wake recedes until levels similai' to those existing 
at the start ol the traverse are established. This behavior 
is directly related to the changes in pressure over the air- 
toil surtace caused by the presence of the wake. As the 
wake moves past the stator, the local tlow undergoes 
transient fluctuations in velocity and turbulence inten- 
sity that affect the boundary layer, and consequently, the 
surface pressure distribution. 

The maximum variation in c, resulting from w-ake-stator 
interaction ctYects can be de lined as 
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In ihc above equation. (Cj )^j, the solid stator c, value at 
the start ol' the wake traverse (t = OX) s). was obtained 
from a converged solution with the wake held stationary 
at a distance of approximately 2 chords above the air- 
foil. is the minimum c, value associated with 

wake convection. Thus, for the solid stator airfoil, 
(Ac,) is29Vf. 

U ma X 



f igure 5 - Idfeci of wake passage on stator lift coeflicient, line 
gnd level. .M = l). 166. Re = 1. 125x1 ()^ V = 78c/s 


Hffect of Surface Porosity on Airfoil Pressure 

[he extent and location of the porous regions will 
depend, of course, on the intended application. For the 
reduction of unsteady interactions in turbomachinery 
environments, the overriding performance concern is the 
minimization of (low separation on the stator surtaces. A 
pai'ametric study involving 45+ different porous conlig- 
urations was performed in order to isolate tho.se that 
would minimize the change in lift cau.sed by w'ake con- 
vection, without compromising airfoil performance. The 
valuable parameters were the extent and location of the 
porous patches (ranging from full-chord porosity to 
small isolated regions near the LH and TH), the conligu- 
ralion of the assumed plenum chamber (single or multi- 
ple. either ctmnected or disconnected), and the surface 
porosity (either constant or tapered with distance along 
the surface). A value of 22'^ porosity^ was chosen as a 
basis. 

The study showed that as long as the porous patch on the 
leew'ard side does not extend past the 107 ( chord loca- 


tion. the flow remains attached over the entire airfoil 
suiiace. It was also lound that, because of airfoil geome- 
try. the porosity on the leeward side had to be tapered 
from the nominal to less than 107 in the region 0.05< x/ 
c < 0.10 in order to preclude the formation of a small 
recirculation area immediately downstream of the 
porous patch. Vaiiations in porosity of ±\07 from the 
basis had a small eflect on airfoil lift. The contiguraiitms 
yielding the largest unsteady lift reductions were those 
that involved communication, through the plenum, 
between regions of high pressure differential ahead of x/ 
c = 0.15. 

Some of the viable configurations obtained during the 
parametric study are presented next. The .series, shown 
in Figure 6, has a leeward side porosity distribution of 
227r in the region 0 < x/c < 0.05, then tapered ellipti- 
cally from 227 to 107 in the region 0.05 < x/c < 0.10. It 
also includes an isolated 227 porosity patch in the 
region 0.95 < x/c < 1 .0 aimed at alleviating TF How sep- 
aration. It has been determined that this patch is not nec- 
essary, since, in general, the viable conligurations do not 
induce (low sepai'ation. However, it was included early 
in the study and kept for some of the configurations. Its 
effect is negligible. The porosity distribution on the 
windward side consists of a 227 porous patch with an 
extension that varies from 107 of the chord up to 257, 
507 and full-chord, for configurations A. B. C. and D. 
respectively. The main porous region in each configura- 
tion is continuous, i.e., it is defined as a single patch 
with a single implicit plenum chamber. 


Confiauration A 


C.'onliizuralion H 


Contiguration C 


Configuration 1) 


Figure 6.- Sample poixius disiribuiions. 
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The effecl of porosity on the How in the vicinity of the 
stator I.E is depicted in Figure 7, for t = 0.0 s. Wake-sta- 
tor interaction effects are negligible. The small amount 
of mass transfer characteristic ol passive porosity is rep- 
resented by velocity vectors whose length* amplified 
several times for visualization purposes, is proportional 
to the magnitude of the tlow normal to the surtace. Cor- 
responding surface pressure coellicient distributions are 
also presented in the ligures. The How ai'ound the solid 
stator TH (Figure 7a) is provided as a reference. 

Porosity alters the effective aerodynamic shape of a 
body by modifying the pressure distribution over its sur- 
face. As expected, these changes affect the How in the 
vicinity of the body. The expected effect of porosity on 
the local How in the vicinity oi the leading edge ol an 
airfoil is the creation of stable regions of inflow (ppienum 
< Pdomam) surfaco prcssuie, and outtlow 

(Ppionun. > Pdc„„a,n) io uTcas of low suftace pressure. These 
high and low pressure areas are usually located on the 
windwai'd and leeward surfaces of the airtoil, respec- 
tively. In regions of inllow, some ol the decelerated tluid 
particles in the boundary layer ai'c removed, resulting in 
a thinner boundary layer. In regions of outllow, addi- 
tional energy is supplied to the tluid particles being 
retarded, resulting in the creation ot a sub-layer near the 
wall that reduces shear forces and displaces the bound- 
ary layer away from the surface. However, the numerical 
results presented here indicate that porosity tends to 
induce a time-dependent oscillatory pattern ol small 
intlow-outllow regions neai' the stator leading edge, 
which is well established be lore wake ellects come into 
play. 


(a) Solid stator 



The causes of the oscillatory pattern are not well under- 
stood. but are most likely related to airftiil geometry and 
(low incidence. As mentioned earlier, an analytical deli- 
nition for the airfoil was not available. Thus, a circular 
arc was fitted to the given discrete points to form a Lh. 
Although the lilting process was carefully done, minute 
slope discontinuities could not be avoided. Because the 
airfoil is subject to off-design conditions, alypically low 
levels of lift are generated at the I.H, as seen in the Cp 
distribution for the solid airfoil (Figure 7a). It is possible 
that the combined effecl of these two lac tors triggered 
the oscillatt)ry inflow-outllow pattern shown here. 

rhe periodic, oscillatory pattern starts at the FH. where 
air always Hows into the plenum, and travels down- 
sti'eam on both leewai'd and windwai’d sides of the air- 
foil. On the leeward side, the pattern travels to the all 
end of the patch and beyond. This additional distance 
depends on the extension ot the patch on the low’er sur- 
face, but is between an exU’a ().()5c for configuration A to 
about an extra 0.20c for conliguration 1). The amplitude 
of the disturbances peaks at the alt end ot the patch (x/c 
= 0.1) for all configurations, and diminishes from there. 
On the windward side ot the airtoil. the disturbances 
travel to approximately ."SO^^ of the patch extension. 
Their amplitude peaks shortly aft of the I.F. and dimin- 
ishes as the disturbances travel downstream. One period 
of this oscillatory pattern is presented in Figure 8, for 
configuration A, The period considered in the figure is 
early during the wake traverse in order to exclude wake- 
stator interaction effects. 



7.- Idfect of surface porosity on Mach number contours and surface pressure distributions lor tlow 
amund the stat()r \M at 1 = 0.0 s. Medium grid, M« = 0, 166, a = 6"^. 
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(e) Configuration I) 



000 005 0 10 


x/c 


1 5 r- 




-0.5 h 



solid 

porous D 


1 5 


0 


0 ^ 

X/C 


i— 

0.75 


1 


I'igure 7.- Concluded. 


Regardless of the unsteadiness caused by the periodic 
inllow-outflow oscillations, it can be seen Irom the pres- 
sure distributions that the How remains attached over the 
entire surface for all configurations except D, for which 
the How over the leewai'd side has slightly separated in 
the TH region. 

The periodic, oscillatory behavior described here would 
probably not arise if the airfoil were placed within a sta- 
tor cascade, and subject to How conditions typical of tur- 
bomachinery environments. However, it has been 
observed experimentally during an internal NASA 
LaRC wind tunnel test of a spherical nosc/cylinder fore- 
body model with equivalent porosity levels-f 


-i.5r 



t - 0.00237 s 

' t - 0,00282 s 

1 - 0.00326 s 

t = 000370s 
t = 0,00423 s 

^ ^ 0 0.25 0.5 0,75 1 

x/c 

f igure 8.- Hffccis of pomsity-induced How unsteadiness on 
surface pressure distribution, configuration A. Medium grid 
level. = 0. 1 66. Re = U 25x I O''. 


KtTcct of M(wing Wake on P(m)us Airfoil Behavior 

It has been observed from surface pressure animations 
(not show'n) that wake-stator interaction tends to reduce 
the amplitude of the porosity-induced oscillations on the 
leeward side of the airfoil. For configuration A, which 
has the least amount of porosity, the strength of the 
interaction is sufficient to cancel the oscillations alter 
the wake intersects the stator LH. The cancelling etteci 
of the interaction decreases as the extension of the patch 
on the windward side increases. For configuration D. 
wake impingement momentarily halts the oscillator) 
behavior over the patch surface. However, it stalls 
immediately aft of the patch, and travels downstream on 
the displaced boundary layer an additional 0.1c. After a 
while, the behavior observed prior to wake approach 
reestablishes itself. 

Wake-stator interaction tends to amplify the oscillations 
on the windward side of the airfoil and to momentarily 
extend their uavelling distance. For configuration A, the 
amplitude of the oscillations is seen to increase and then 
decrease as the wake intersects and convects past the 
stator LH. The amplitude remains slightly higher than 
before wake approach, but decreases as the waves travel 
downstream. The distance travelled by the disturbances 
extends to about x/c = 0. 15 (0.05c beyond the end of the 
patch) immediately after the wake intersects the stator 
LH, eventually shrinking back to hall the patch length. A 
similar behavior is observed for the other configura 
tions, although wave amplilication seems to be propor- 
tional to the extension of the patch on the windward 
side. For configuration D, the oscillations are consider- 
ably amplified immediately after closest approach, 
momentarily travelling downstream over the entire 
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length ot the patch. The amplitude diminishes somewhat 
with time, and the travelled distance is reduced to halt 
the patch length {0.5 x/c). 


Hllect of Su rface Porosity on Airtnil Unsteady 1 iO 


The ellect ot wake passage on porous airfoil C| is 
depicted in Figure 9. Note that, as lor the solid case, the 
c, is seen to decrease and then increase as the wake 
approaches and recedes from the stator, respectively. 
Several chajacleri sties can be iden tilled in the time his- 
tories, all direct consequences of the changes in effec- 
tive geometry and increments in local Mow unsteadiness 
brought about by porosity. 
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I iguic 9, Fitl ci)clticient time histories for various porous air- 
toil cimtigurations, medium grid level. M„ = (I 166. Re = 
l.l2,‘>xll»f>.V = 78c/,s 


Note Irom Figure 9 that the converged lift coeflicient at 
the start of the traverse (t = 0.0 s) is different for all con- 
ligurations. This is to be expected, since it is ass(x:iated 
with different effective aerodynamic shapes. The aver- 
age deciement in initial C| ranges Irom approximately 
4.57r to about 107 for contigurations D and C, respec- 
tively. The average increment in airfoil drag coefticient 
at the stall ot the traverse (not shown) ranges from about 
17 for conliguration A to a drastic 507 Vor ctmligura- 
tum D. The unsteadiness inherent in the tlow as a result 
ol porosity, which is represented by the amplitude of the 
small oscillations present in the time histories, appeal's 
to be pioportional to the extent ol the porous reeion. 
This trend, in conjunction with the alterations in effec- 
tive geomeliy, is believed to be the cause of the changes 
in airtoil lilt and drag previously noted. 


It can also be observed from Figure 9 that, for the global 
period under consideration ('-0.05 s), the lift coeflicients 
associated with the porous contigurations do not recover 
to the levels existing at the start of the wake traverse, 
unlike those lor the solid airfoil. In general, this hystere- 
sis ellect seems to be proportional to the extent of poros- 
ity over the lower surface, and it may eventually 
disappear il longer periods are considered. 

An attempt was made during the study to compare 
porosity effects on the basis of equal lift. To do so, the 
angle ol attack of conliguration A was increased until its 
Cj value at the start of the traverse matched the value cor- 
responding to the solid airtoil. Since the increment in 
angle ol attack was relatively large (approximately 2.5*9. 
the resulting increment in drag was prohibitive, and the 
ability of porosity to reduce unsteady lift was dimin- 
ished. It is apparent that, in order to compare on the 
basis ol equal lilt/drag, the shape ol the airfoil would 
have to be redesigned/optimi/ed including porosity. 
Such an undertaking is outside the scope of the present 
investigation. 

Regardless of the changes in effective geometry and 
local Ilow', all the porous configurations being consid- 
ered were able to reduce (Acj . This reduction, ref- 
erenced to ranged Trom 167 to 267 for 

configurations B and D, respectively. Considering the 
changes in lift and drag coefficients at t = 0.0 s, (Ac,) 
and (Ac\,)j^ . respectively, it appears that configuration A. 
with (Ac,)^, - -57 and (Ac;,)^^ ^ 17. is the best per- 
lormer with a reduction in unsteady lift of approxi- 
mately 187 for the medium grid level results shown 
heie. A fine grid solution, obtained lor configuration A 
only, indicated that the reduction in unsteady lift could 
increase up to 217 with similai' performance losses. 


Concluding Remarks 

The potential ot passive porosity technology as a mech- 
anism to leduce the unsteady lilt on airfoils exposed to 
wake convection has been evaluated. To do so. a typical 
tan stator airfoil was immersed in a free field and sub- 
jected to the ellects of a transversely moving wake. 
Solutions were obtained for a solid airfoil bascHne, and 
lor a series of porous configurations in order to identify 
those that would reduce unsteady lift without compro- 
mi.sing airfoil performance. It was found that as long as 
the extension of the porous region on the leewai'd side 
did not exceed 107 ol the chord, the tlow over the entire 
airfoil surface remained attached. Despite a porosity- 
induced unsteady pattern of small inflow-ouillow 
regions neai' the stator leading edge, it was determined 
that the use of porosity could reduce the magnitude of 


11 

American Institute of Aeronautics and Astronautics 



the unsteady lilt resulting I'rom wake-stator interaelion 
upwards of 187^, without incurring significant perfor- 
mance losses. 

The work presented here constitutes the first stage ot a 
study designed to assess the benefits of using passive 
porosity to reduce interaction noise in turbomachinery 
environments. The second stage ot the study, which 
quantities the reductions in unsteady lift in terms ot 
changes in radiated rotor-stator interaction noise, will be 
presented in a future paper. 
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